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We use polarized magneto- Raman scattering to study purely electronic excitations and the 
electron-phonon coupling in bulk graphite. At a temperature of 4.2 K and in magnetic fields up 
to 28 T we observe ii^-point electronic excitations involving Landau bands with A|n| = and with 
A|n| = ±2 that can be selected by controlling the angular momentum of the excitation laser and of 
the scattered light. The magneto-phonon effect involving the E2g optical phonon and iT-point inter 
Landau bands electronic excitations with A|n| = ±1 is revealed and analyzed within a model taking 
into account the full kz dispersion. These polarization resolved results are explained in the frame of 
the Slonczewski-Weiss-McClure (SWM) model which directly allows to quantify the electron-hole 
asymmetry. 

PACS numbers: 73.22.Lp, 63.20.Kd, 78.30.Na, 78.67.-n 



I. INTRODUCTION 

Magneto-optical spectroscopy has been used to study 
the graphene/graphite structures for long time, but 
mostly transmission/reflectivity type of experiments (es- 
sentially in far-infrared spectral range) have been ex- 
plored so far— Relevant for these systems, low en- 
ergy electronic excitations can be optionally probed with 
Raman scattering methods^^. These methods have been 
widely used to investigate phonon resonances in different 
carbon materials but advantages of combining them 
with application of magnetic fields have been recognized 
only recently—"—. 

Indeed, magneto-Raman scattering experiments pro- 
vide a spectacular demonstration of hybridization be- 
tween optical phonon and electronic excitations in epi- 
taxial graphene^^ and graphene locations on graphite 
surfac Q -^ T^^ . What is perhaps even more important is 
that such experiments allow also to probe purely elec- 
tronic excitations in these systems^^. When a magnetic 
field is applied across the layer, the electronic bands con- 
dense into Landau levels whose characteristic energy lad- 
ders (fan charts) reflect the specific dispersion relation 
of electronic states in the absence of the magnetic field. 
Magneto- Raman scattering can be used to trace selected 
inter Landau level excitations. This method, with all its 
advantages of a visible optics technique (spatial focus- 
ing, polarization resolved measurements), appears now 
as a valuable option for Landau level spectroscopy to 
study other sp^-bonded carbon allotropes, such as bulk 
graphite, investigated in this work. 

In spite of recent efforts^^, the magneto- Raman scat- 
tering response of bulk graphite has not been so far 
clearly identified in experiments. Theoretical studies ad- 
dressing this specific problem seem to be also missing. 



Relevant for our results are, however, predictions for a 
bilayer graphene which is the basic unit in construc- 
tion the Bernal-stacking of graphite. 

As a layered material, graphite is strongly anisotropic 
but still represents a three dimensional electronic sys- 
tem. In this respect, graphite is very different from its 
layer-components: purely two-dimensional monolayer or 
bilayer graphene. Electrons can easily move within the 
graphitic layers, but electronic states also remain disper- 
sive in the direction across the layers. When a magnetic 
field is applied across the layers, we deal with Landau 
bands in graphite in contrast to discrete Landau levels 
in purely 2D systems (monolayer or bilayer graphene). 
Band structure of graphite is commonly described using 
the SWM mode l^^i^^ . There is a general consensus on the 
validity of this tight binding approach. However, it im- 
plies up to seven hoping-integral parameters and the ex- 
act values of some of them remain controversial and this 
is despite very long history of the research on graphite. 

Our magneto-Raman scattering experiments on 
graphite reveal a series of features which we attribute to 
electronic excitations between Landau bands in this ma- 
terial. The magnetic field evolution of the phonon 
excitation is also investigated. The observed electronic 
and phonon excitations are shown to follow the charac- 
teristic selection rules defined by the appropriate transfer 
of angular momentum between the incoming and outgo- 
ing, circularly-polarized photons. Polarization-resolved 
experiments together with appropriate modelling of elec- 
tronic bands allow us to conclude about the relevant 
electron-hole asymmetry of the graphitic bands. Al- 
though our experiments are mostly sensitive to electronic 
states in the vicinity of the particular K-point of the 
graphite Brillouin zone, the three dimensional nature of 
this material is clearly reflected in the spectral shape of 
the observed electronic excitations and in the character 
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of the investigated magneto-phonon resonances. Notably 
the phonon of graphite is shown to hybridize with a 
quasi-continuous spectrum of inter Landau band excita- 
tions, in contrast to similar effect in graphene involving 
the discrete LL transitions. 

The paper is organized as follows: section II describes 
the experimental set-up, while the electronic properties 
of bulk graphite and the Raman scattering selection rules 
are presented in section III. Experimental results ob- 
tained in the co- and crossed-circular polarization con- 
figuration are presented in section IV, together with the 
magneto phonon effect of bulk graphite and experiments 
performed at room temperature. We finally present our 
conclusions. 



II. EXPERIMENT 

The polarized Raman scattering response of bulk 
graphite at low temperature and in high magnetic fields, 
was measured with a home made confocal micro-Raman 
scattering set-up. A monolithic miniaturized optical ta- 
ble, made of titanium, has been designed to operate in a 
helium exchange gas at T= 4.2K and in high magnetic 
field environments. A mono- mode optical fiber with 5 /am 
core was used to bring the A = 784 nm excitation from 
a TiiSaphire laser to the sample. The excitation beam 
is focused with lenses down to ~ 1 /im spot. Scattered 
light is then collected by a multi-mode 200 /im optical 
fiber before being analyzed by a 500 mm spectrometer 
equipped with a diffraction grating and a nitrogen cooled 
CCD. Our miniaturized optical table can host different 
optical filters that are used at liquid helium temperature 
and in high magnetic fields. Band-pass filters (laser line 
and notch filters) are used first to clean the laser coming 
out of the mono-mode fiber and to filter the elastically 
scattered laser before the collection optical fiber. They 
impose a cut-off energy of ~ 350 cm~^ from the laser 
line. A set of quarter wave plates and linear polarizers, 
placed close to the sample, are then used to circularly 
polarize the excitation beam and the collected signal to 
achieve both co- and crossed-circular polarization con- 
figurations. The two different co-circular and 
cr+/(7+) or cross-circular (g~ j and j g~) configu- 
rations were obtained by changing the direction of the 
magnetic field with respect to the light propagation di- 
rection. Optical power on the sample was fixed to 5 mW. 
The sample is mounted on a set of translation piezzo 
stages which allow to move the sample under the laser 
spot with a sub-micrometer resolution. 

The surface of natural graphite is known to be inho- 
mogeneous and graphene flakes decoupled from the sur- 
face can be found. They have been evidenced by low 
temperature STM experiments^^, studied by an EPR-like 
technique in magnetic fields^^ and mapped with micro- 
Raman scattering measurements in high magnetic fields 
revealing their graphene-like electronic excitation spec- 
trumi^ and the associated magneto-phonon effect i^:^. 



To unravel the electronic properties of bulk graphite, we 
place the laser spot outside of these decoupled graphene 
flakes and measure the magnetic field evolution of the 
Raman scattering spectrum. Similar results were ob- 
tained using two different sources of natural graphite and 
Highly Oriented Pyrolitic Graphite (HOPG) and in the 
following, we will only discuss the case of natural graphite 
which shows, in our experiment, a slightly higher signal 
level. 



III. THEORETICAL OUTLINE 
A. Graphite band structure 

A conventional description of the band structure of 
bulk graphite and its evolution with the magnetic field 
relies on the SWM model with its seven 70,..., 75, A 
tight binding parameters22a^iSP . This model has been 
used to describe most of previous data obtained from 
magneto-transportf^i — ^ infrared magneto-reflectivity^^, 
and magneto-transmission^^ ^^^^^ experiments. It predicts 
the existence of massive carriers near the K point with 
a parabolic in-plane dispersion and of massless carriers 
near the R point with a linear in-plane dispersion. The 
Fermi energy is ~ 15 meV and, under an applied mag- 
netic field. Landau bands are formed with a continuous 
dispersion along /c^, from equally spaced and linear in 
B Landau levels at the K point to non-equally spaced 
and \fB evolving Landau levels at the 11 point^^. Even 
though there is still no consensus concerning the pre- 
cise values of the SWM parameters, mainly because of 
the different energy ranges probed in different experi- 
ments and because of the lack of polarization resolved 
measurements that would reveal unambiguously the ef- 
fect of electron-hole asymmetry, the validity of the SWM 
model is generally accepted. 

The evolution in magnetic field of electronic levels in 
bulk graphite has been calculated following the approach 
of Nakac^^ and Landau levels are labelled following the 
bilayer graphene convention^^ as sketched in Fig. [TJ The 
infinite order magnetic field Hamiltonian was reduced to 
a 600 X 600 matrix before the diagonalization procedure. 
Electronic excitations are labelled as L^^n' where n and 
n' are the indices of the Landau levels involved in the 
excitation n ^ n' . 

Instead of the full SWM model, it is often sufficient 
to use the effective two-parameter model, which gives 
parabolic dispersion in the plane with the slope of the 
parabolas depending on A:^, the wave vector measured in 
the units of the inverse inter-layer spacing {kz = at the 
K point, and kz = 7t/2 at the H point). This model is 
obtained by (i) neglecting all SWM parameters except 
two, 70 and 71, the intralayer and the inter layer nearest- 
neighbor hopping integrals, respectively, and (ii) project- 
ing the resulting /c^-dependent 4x4 Hamiltonian on the 
two low-energy bands. At each value of /c^, the Hamilto- 
nian is identical to that of a graphene bilayer determined 
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FIG. 1. Evolution of the Landau levels at the K point of 
the graphite first Brillouin zone calculated within the SWM 
model and Landau levels' indices. 



by the effective parameters 70 and 7^ = 271 cos kz . Since 
at the K point /c^ = 0, the corresponding 7^ is twice en- 
hanced with respect to 71 describing the real graphene bi- 
layer^^'^^. This effective two-parameter parabolic model 
has proven to bring a fair frame to describe magneto- 
optical experimental^. 

Still, the effective two-parameter parabolic model 
misses several effects, such as the in-plane trigonal warp- 
ing, described by the 73 parameter of the full SWM 
model, as well as the electron-hole asymmetry. The latter 
is contributed by all remaining parameters (72, 74, 75, A). 
To simplify the analysis, one can notice that at kz =0, 
the 72, 75, A do not enter independently; eigenstates and 
transition energies depend only on their combination 
Aeff = A+275-272, as can be seen from Eqs. (jA3)) . (|A4|) 
of the Appendix. This combination determines the asym- 
metry in the positions of the split-off bands exactly at 
the K point. For non-zero values of the in-plane wave 
vector, the electron- hole asymmetry is determined by 
both 74 and Ag//, and it is not easy to separate the 
effects of the two. Throughout the paper, we use a set of 
SWM parameters derived from our polarization-resolved 
magneto- Raman scattering experiments. We will detail 
in the following sections how these parameters are de- 
termined, essentially from the crossed circular polariza- 
tion configuration measurements which allow to directly 
quantify the electron-hole asymmetry. 



B. Selection rules for Raman scattering 

Generally, the Raman scattering selection rules for 
graphite can be deduced from its symmetry group, Dq^- 
However, the low-energy electronic Hamiltonian, as de- 
rived from the SWM model, is dominated by terms which 
have a higher symmetry. As a consequence, the electronic 
dispersion is almost isotropic with respect to continuous 
rotations of electronic momentum around the H — K — H 
line. This isotropy is broken only by the trigonal warp- 
ing term, governed by the 73 parameter of the SWM 



model. This term restricts the symmetry to the momen- 
tum rotations by ±27r/3. As a result, the quantum num- 
ber m^, associated with the rotational symmetry (angu- 
lar momentum), which could assume all integer values 
from —00 to 00 in the case of the continuous rotational 
symmetry, has only three distinct values for the three- 
fold symmetry: all m^'s differing by a multiple of 3, be- 
come equivalent. However, since the trigonal warping is 
small at energies, relevant for the present work, all tran- 
sitions which are allowed by Dq^, can be further classified 
into strongly allowed (those which are allowed even when 
73 = 0, thus surviving the continuous rotational symme- 
try), and weakly allowed (those which require a non-zero 
trigonal warping). The wording "strongly" and "weakly" 
allowed does not, however, refer to the apparent intensity 
of the excitation. 



An example of a weakly allowed process is the Ra- 
man process responsible for the G band feature. Indeed, 
trigonal warping is necessary for the appearance of the 
G band feature, as was shown for monolayer graphene 
in Refs. [s^ andlssi: in bilayer graphene and graphite the 
situation is analogous. The reason is that the circularly 
polarized photons carry an angular momentum rriz = ±1, 
and out of two degenerate phonon modes, one can also 
choose two linear combinations carrying angular momen- 
tum ruz = ±1. Then if a cr+ photon (m^ = —1) is ab- 
sorbed by the sample, the only way to conserve angular 
momentum is to emit a a~ photon {rriz = +1) and a 
phonon with = +1. Then the total change in the an- 
gular momentum is +3, which is equivalent to 0, thanks 
to trigonal warping. Thus, the G band feature is observed 
in the cross-circular polarization configuration. 



An inter-Landau-level excitation n~ (n + m^)^ car- 
ries an angular momentum m^. This immediately desig- 
nates the transitions with A|n| = and A|n| = ±2 as the 
strongly allowed processes in the co-circular and cross- 
circular polarization configurations, respectively. As in 
the case of a graphene monolayer, the strongly allowed 
coupling for the magneto-phonon effect in a graphene 
bilayer^^ involves the optical phonon and electronic exci- 
tations with A|n| = ±1 (optical-like excitations). These 
optical-like excitations are not expected to be Raman 
active without trigonal warping. Nevertheless, if 73 is 
included then the trigonal warping is expected to weakly 
allow these excitations. For instance, an excitation 
Ln-^{n-\-i)+ with A|n| = +1 can be observed in Raman 
scattering experiments, because of the trigonal warping 
induced mixing of the levels n~ and (n + 3)~, through 
the excitation L(^^3)- which is strongly allowed in 
Raman scattering. Interesting is the fact that A|n| = ±1 
are expected to be observed in the same crossed circu- 
lar polarization as A|n| = ^2 excitations. A signature 
of such transitions has indeed been seen in the present 
experiment. 



IV. EXPERIMENTAL RESULTS AI 
DISCUSSION 

A. Co-circular configuration 

Fig. [2fa) shows representative Raman scattei 
tra in the co-circular polarization configurat: 
the B = spectrum subtracted. In this confi 
the G band feature is absent, as discussed in ^ 
Starting from 5 ~ 2 T, many magnetic field d 
features are observed, with energies increasing 
magnetic field. The observed line shape is stron 
metric with a long tail at high energies due 
that the states involved in the scattering proce 
belong to discrete Landau levels, but to Land 
with a dispersion along kz. The width of these 
about 50cm~^, is mostly determined by this / 
sion, rather by any scattering mechanism. To ^ 
quantitative understanding of the line shape, 
that the Raman matrix element does not depei 
energy of the electronic states involved^^, so the 
intensity I{uj) is simply proportional to the j 
sity of states for excitations between Landau 
and 
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where e^± are the Landau level energies, and 
netic length, Ib = (e5/c)~^/^, keeps track of th 
level degeneracy. In the two-parameter parabo] 
the Landau level energies are given by 
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where a is the in-plane distance between neighboring car- 
bon atoms. For e^± given by Eq. ([2]), Eq. ([T]) gives 
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FIG. 2. a) Representative Raman scattering spectra with 
subtracted B = T spectrum in the a~ /a~ configuration 
showing many asymmetric magnetic field dependent features. 
Spectra are shifted for clarity. The gray region indicates an 
energy range where a magnetic field dependent background 
feature affects the normalized spectra, b) Raw Raman scat- 
tering spectra to compare, at 3 different values of the magnetic 
field, the scattered intensity. Solid and dashed lines are cal- 
culated within the two-parameter parabolic model (see main 
text for details). The vertical dotted line represents the cut- 
off energy of the optical edge-filter used to remove the laser 
line. 



where uOn = — e^-|^ ^q. This expression describes 
well the tail on the high-frequency side of each peak, and 
has a square-root-type singularity at ^ cj^. The singu- 
larity can be cut off by replacing the (5- function in Eq. ^ 
by a Lorentzian with full width at half- maximum F. This 
results in a somewhat more complicated, but still explicit 
analytical expression for the Raman spectrum: 
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In fact, the simple expression ([2]) does not reproduce well 
the peak positions; the full SWM model taking into ac- 
count the electron- hole asymmetry is needed for this, as 



will be discussed below. Still, if we take Eqs. (j4]), ^ and 
substitute for Un the actual peak positions, the result- 
ing expression describes the spectrum remarkably well, 
as shown in Fig. 2(b) for two values of the magnetic 
field. We have taken F = 14 cm"^ for B = 12T and 
F = 20cm~^ for B = 20 T. The discrepancy between 
Eq. (jll and the experimental spectrum around the min- 
ima between the peaks can be improved by taking the 
fun SWM result for e^±. 

Fig. Efb) shows that there exist an electronic contri- 
bution to the B = Raman scattering spectrum of bulk 
graphite, which is flat as a function of energy (between 
300 and 1200 cm~^), but that can be identifled by apply- 
ing a magnetic fleld. With the three Raman scattering 
spectra presented in this flgure, one can directly compare 
the amplitude of the scattered light. When a magnetic 



TABLE I. SWM parameters used for the calculations 



SWM parameter 


full model (eV) 


parabolic model (eV) 


70 


3.08 


3.15 


71 


0.38 


0.38 4) 


73 


0.315 





74 


0.044 





A + 275 - 272 


0.22 






field is applied, the energy independent response from 
low energy electronic excitations transforms into discrete 
features due to Landau quantization and one can iden- 
tify the " apparent background" signal in this experiment 
with the scattered amplitude at energies lower than the 
1,2-^2+ excitation, which is set to zero in Fig.[2fb). In this 
particular experiment, we use magnetic field to change 
the electronic excitation spectrum, but one could expect 
a similar trend for gated graphene flakes. Changing the 
gate voltage would result in a continuous change of the 
2Ef threshold in analogy to the tuning of the absorption 
threshold as a function of the Fermi energy observed in 
infrared absorption experiments^^. One may also notice 
that the amplitude of the scattered light between any 



two given peaks L^ 



and L 



(n+l)-,(n+l)H 



does not de- 



pend on the magnetic field. This can be seen directly 
from Eq. ([!]), which does not rely on approximation (j2j): 
as long as Landau level energies are proportional to B 
(which is the case for uj < 1000 cm~^), the magnetic field 
drops out of /^(cJn+i)- 

The upper panel of Fig. [3] shows a false color map of 
the scattered intensity as a function of the magnetic field 
for the two a~ /a~ and (j+/cr+ configurations. The evo- 
lution of these features is nearly linear with the mag- 
netic field which is expected for carriers with a parabolic 
dispersion. The evolution of the maxima of the scat- 
tered light intensity as a function of the magnetic field 
is presented in the lower panel of Fig. [3] for the a~ /a~ 
(open dots) and j (black dots) configurations. We 
present in this figure the results of the corresponding ex- 
citations calculated within the two-parameter parabolic 
model with 70 = 3.15 eV, 7^ = 0.760 eV (black lines). In 
a first approximation, this model allows to describe well 
the overall behavior of excitations involving Landau lev- 
els with n > 2. The clear limitations of this model appear 
at low energies with the 1/2-^2+ excitation which is over- 
estimated, mainly because of the missing 73 parameter 
relevant at low energies. 



B. Crossed circular configuration 

The crossed circular polarization configuration is ex- 
pected to select excitations of the E2g symmetry which 
includes the optical phonon (G band feature) and elec- 
tronic excitations involving hole and electron states with 
A|n| = ±2 (Ref. [l7| ), which correspond to angular mo- 
mentum rriz = ±2. We present in the upper panel of 
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FIG. 3. Upper panel: False color map of the scattered inten- 
sity at T = 4.2 K as a function of the magnetic field in the two 
co-circular configurations (raw data). Lower panel: Position 
of the maxima of scattered light as a function of the magnetic 
field in the /a^ (black dots) and (t~ /(t~ (open dots) con- 
figurations. Solid black lines are the A|n| = electronic exci- 
tations calculated within the two-parameter parabolic model 
using the parameters described in Table. I. 



Fig. m a false-color map of the scattered intensity as a 
function of the magnetic field for both crossed circular 
polarization configuration and, in the lower panel of the 
same figure, a false color map of the derivative with re- 
spect to the magnetic field (A5 = 0.32 T and spectra 
acquired every 0.08 T) of the same data. In these fig- 
ures, many magnetic field dependent features can be ob- 
served, as well as a pronounced oscillatory behavior of 
the G-band feature that will be discussed in detail in 
the next section. All of these features have an almost 
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FIG. 4. Upper panel: False color map of the scattered in- 
tensity at T=4.2K as a function of the magnetic field in the 
two crossed circular configurations. Lower panel: Derivative 
of the scattered light intensity as a function of the magnetic 
field. 



linear dependence on the magnetic field confirming that 
i^-point carriers are involved. 

Most of the excitations shown in Fig. |4] are expected 
to follow the A|n| = ±2 selection rule. What is surpris- 
ing, especially in the j g~ configuration, is the seeming 
coincidence of these excitations with the A|n| = ±1 ex- 
citations which we attribute to the effect of the trigonal 
warping discussed in Sec. IIIIBI and which couple to the 
E2g phonon. Indeed, a closer inspection of Fig. |4] shows 
that the magnetic field dependent features that are ob- 
served are not those that couple to the phonon. This 
is particularly visible in the cr~ / cr^ configuration where 
the electron-phonon interaction is resonant at 5 = 21.8 T 
while the electronic excitation crosses the phonon energy 
at B = 24.1 T. The coincidence of the magnetic field 
values corresponding to resonant electron-phonon inter- 
action revealed by the oscillations of the phonon and the 
magnetic field values at which the visible electronic exci- 



tations cross the phonon energy is a result of the electron- 
hole asymmetry. This asymmetry leads to nearly degen- 
erated (n + l)~ (n — 1)+ Raman active excitations and 
(n — 1)~ optical-like excitations. This effect is less 

pronounced for lower values of the Landau level index 
which makes this particular resonance at high fields well 
separated from the Raman active electronic excitation. 

In the following, we search for appropriate parame- 
ters of the SWM model, taking into account that (i) 
A|n| = ±1 electronic excitations should cross the phonon 
energy at the anti-crossing points and (ii) that the ob- 
served electronic excitation should correspond to the 
A|n| = ±2 electronic excitations in both crossed circu- 
lar polarization configurations. These two requirements 
are indeed fulfilled by the parameters summarized in Ta- 
ble I. The results for different electronic excitations de- 
termined from the effective SWM model together with 
the observed excitations in both crossed circular polar- 
ization configurations are presented in Fig. [5l In this 
figure, the gray and black lines correspond to the calcu- 
lated excitations in the g~ j and j g~ configurations 
respectively. The open triangles and open squares cor- 
respond to the maxima of scattered intensity observed 
in and <j+/<j~ configurations respectively. The 

gray and black circles correspond to resonant magnetic 
fields for the magneto-phonon effect in these two different 
configurations. As can be seen in Fig. O the set of SWM 
parameters presented in Table I describes well most of the 
observed excitations both in the co-circular configuration 
for which the electron-hole asymmetry does not influence 
the different excitations energy and in the crossed circu- 
lar configuration where a strong electron-hole asymmetry 
is observed. 

To determine these parameters, we start from the 70 
and 7i values determined from magneto-transmission ex- 
periments on similar samples and from the 73 and 
74 values proposed in Ref. [sl. We then introduce 
the electron- hole asymmetry by increasing Ag// starting 
from values presented in Ref. 38 and compare the results 
with experimental data. The value of Ag// is then grad- 
ually increased in order to describe the observed energy 
difference in both crossed circular configurations. Intro- 
ducing the electron-hole asymmetry has an effect on the 
70 value which have to be slightly decreased (by 4%) to 
describe both crossed and co-circular polarization exper- 
iments while 7i, 73 and 74 are kept constant. It should 
be noted that, even though we have reduced the num- 
ber of parameters by introducing Ag//, there is not a 
unique pair of parameters 74 and Ag// able to describe 
our data as both these parameters affect the asymmetry. 
The one we present in Table I keeps 70 within the ac- 
cepted values and allows comparison with other works. 
For the SWM parameters proposed by Brandt et al?'^ ^ 
^eff = 0.108 eV, which is, for the same 74 value, twice 
lower than the Agjj needed to reproduce experimental 
results in the present work. 

The parameters presented in Table I are sufficient to 
describe i^-point carriers but are not valid for /c^ 7^ 
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FIG. 5. Maxima of scattered intensity as a function of the magnetic field for a~ /a^ (gray open triangles) and /(t~ (black 
open squares), together with A|n| = =b2 electronic excitations calculated using the parameters presented in Table. I within the 
two-parameter parabolic model (black and gray solid lines) and A|n ss). 
The gray and black circles correspond to the phonon resonances o ely. 
The line L2- is calculated at /c^^O.S. 



where the seven SWM parameters need to be specified. 
One possible combination is to set A and 72 to their stan- 
dard values A = -0.008 eV and 72 = -0.02 eV^^. This 
leads to 75 = 0.094 eV and of course, at the i^-point, 
this set of parameters brings exactly the same excitation 
spectrum as the one calculated using Ag// = 0.22. These 
parameters can be used to calculate Landau bands ener- 
gies at kz 7^ 0. 

Most of the observed excitation fulfill the Raman scat- 
tering selection rules presented in the previous sections, 
but we also observe traces of optical-like excitations with 
A|n| = ±1: Li,2+ in the g~ j configuration and i^3-,2+ 
in the j g~ . 

All electronic excitations observed in the crossed po- 
larization configuration are significantly weaker than the 
symmetric lines observed in co-circular configuration, but 
unexpectedly, in the g~ j polarization configuration 
two particular features involving the n = 0, 1 Landau 
levels are of much stronger intensity. These two features 
are presented in Fig. [6] (gray curves). They have a sig- 
nificantly different line shape, I/i,3+ being asymmetric as 
the A|n| = excitation discussed in the previous sec- 
tion, while the ^0,2+ feature is much broader. Probably, 
the difference between these two features and the rest of 
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FIG. 6. a) Representative Raman scattering spectra at differ- 
ent values of the magnetic field in (t~ / (black curves) and 
/a~ (gray curves) configurations showing different mag- 
netic field dependent features. 

(n — 1)~ (n + l)+ transitions can be explained by the 
different kz dispersion of the n = 0, 1 Landau levels. 

In the cr^ / cr~ polarization configuration A|n| = —2 
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FIG. 7. Schematics of the observed electronic excitations in 
the two crossed circular configurations and in co-circular con- 
figuration in natural graphite. Observed optical-like excita- 
tions are indicated by dashed arrows. 

and possibly A|n| = +1 excitations are expected. The 
feature with the lowest energy in this configuration is at- 
tributed to Li^2+ excitation at = 0. The first feature 
of the A|n| = —2 series is the L2-,o excitation which 
should be blocked at /c^ = due to the full occupation of 
the n = level, but /c^ 7^ excitations can be observed. 
We attribute the broad feature with a much reduced am- 
plitude as compared to the others (see Fig[6]) to the L2-,o 
excitation at kz = kp ^ 0.3 (dotted lines in Fig. [5j). The 
second excitation with the same change of Landau band 
modulus that could be expected is i^3-,i, which should 
also be blocked because of occupation effects in the n = 1 
level. This feature is not clearly observed in the present 
experiment. 

To summarize all the obtained results, we present in 
Fig. [7] a schematic picture of all the observed excitations 
in the different crossed- and co-circular polarization con- 
figurations. In the cr~ / cr^ configuration electronic exci- 
tations with A|n| = +2 are observed and traces of the 
optical-excitation L3- 2+. The two excitations involv- 
ing the 0, 1 Landau levels also fulfill this index rule. In 
the cr+/(7~, A|n| = —2 excitations are observed together 
with the optical-like excitation 2+ and of the symmet- 
ric excitation L2-,2+- The A|n| = ±1 optical-like elec- 
tronic excitations series are not clearly observed in this 
configuration but the coupling of these excitations with 
the E2g phonon has a strong influence on the phonon 
feature. 



C. Magneto-phonon effect in bulk graphite 

To our knowledge, the magneto-phonon effect involv- 
ing K-point massive carriers in bulk graphite and the 



E2g phonon has never been explored theoretically. This 
situation is of particular interest because of the nature 
of the electronic states involved in the coupling. The 
magneto-phonon effect in graphenei^^^i^ or in bilayer 
graphene^^ involves discrete electronic states which re- 
sults in a fully resonant coupling with the phonon. In the 
case of iC-point carriers in graphite, because of the 3D 
nature of bulk graphite and of the associated dispersion 
along kz^ the electronic states belong to Landau bands 
and the electronic excitation is spread over a wide range 
of energy (line width ~ 50 cm~^ in the case of symmetric 
lines in co-circular configuration) and the interaction is 
also spread on such a wide range of energy. As a result, 
we expect an effect less pronounced than in graphene and 
with asymmetric evolution of the parameters describing 
the phonon feature. 

As already noticed in the previous section, the E2g 
phonon feature shows pronounced oscillations and even 
an anti-crossing behavior when optical electronic excita- 
tions A|n| = ±1 (hardly seen in the Raman scattering 
experiment) are tuned in resonance with the phonon en- 
ergy. This can be seen clearly in the lower panel of Fig. |4j 
We present in Fig. [8] raw Raman scattering spectra in the 
optical phonon range of energy for selected values of the 
magnetic field. The phonon feature energy and FWHM 
evolves a lot with magnetic field and at 5 = 28 T a blue 
shift is clearly visible together with a strong broadening 
of the feature. It is possible to use Lorentzian curves 
to described these spectra, as shown by the dashed gray 
curves in Fig. [HI The parameters (namely, the peak po- 
sition and its full width at half maximum) deduced from 
a Lorentzian fitting of the phonon feature as a function 
of the magnetic field are presented in Fig. [9l^a) and (b) 
for the two crossed polarized configurations. The oscil- 
lations can be observed for magnetic fields higher than 
5 ~ 5 T. One can note on Fig.[9l^a) and (b) that the min- 
ima in FWHM or in Raman shift of the phonon feature 
in the two opposite crossed circular polarizations do not 
occur at the same values of the magnetic field. This dif- 
ference is due to the electron-hole asymmetry in graphite 
discussed and determined in the previous section. 

To describe the magneto-phonon effect theoretically, 
we adopt the same approach that was successfully 
used for the magneto-phonon effect in the monolayer 
graphenei^i2£. Namely, the shift and the broadening of 
the phonon are given by the real part and twice the imag- 
inary part of the complex root uo of the equation 

uj''-ujI = 2uoX [Ub (o; + iV) - Ub=o {^0 + ^r)] , (6) 

where IIb{oj) is the polarization operator, which depends 
on the magnetic field and frequency. We have subtracted 
its value at 5 = 0, so the shift is measured with re- 
spect to cjo, the phonon frequency at zero field (about 
1583 cm-i). 

The simplest two-parameter parabolic model, used in 
Sec. IIV Al to describe Raman scattering due to A|n| = 
electronic excitations, is not sufficient to describe the 
magneto-phonon effect. Four-band approximation is 
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FIG. 8. Raman scattering spectra at selected values of the 
magnetic field in the a~ /a^ polarization configuration (solid 
black curves) and lorentizan curves (dashed gray curves). 
Curves are shifted for clarity. The vertical dashed line in- 
dicates the B=0 position. 



needed, and transitions between levels involving the two 
split-off bands also need to be taken into account. The 
only approximation which is made here is to neglect the 
trigonal warping by setting 73 = 0. For each value of 
kz, and of the Landau level index n, it is then suffi- 
cient to diagonalize a 4 x 4 matrix, instead of an infinite- 
dimensional one, which would be necessary for 73 7^ 0. 
This significantly improves the computation efficiency, 
but introduces a certain error, that will be discussed 
in the following. Relegating the explicit expressions to 
Appendix A, we plot Keuj and 2lmuj for the solutions 
of Eq. dnj by solid lines in Fig. [9l We have adjusted 
A = 3.2 X 10~^ and F = 44cm~^ to reproduce the ampli- 
tude of the oscillations and the degree of asymmetry. We 
added a constant broadening of 5 cm~^ to 2Imcc; which 
corresponds to a broadening of the phonon by scattering 
mechanisms other than the electron-phonon interaction. 
We have used the SWM parameters listed in the left col- 
umn of Table I, and set 73 = 0. 

Setting 73 = changes the inter-Landau-level excita- 
tion frequencies, so their resonance with the phonon oc- 
curs at slightly different values of the magnetic field and 
this explains why the minima and maxima of the theo- 
retical and experimental curves are horizontally shifted 
in Fig. [9l In principle, one could remove this discrepancy 
by slightly readjusting other SWM parameters. However, 
this would introduce a correction to the transition matrix 
elements, and it is not clear whether they would be closer 
to their true values. Thus, we prefer to stick to the pa- 
rameters from Table I, and consider the discrepancy be- 
tween the resonance positions as an order-of-magnitude 
estimate of the error in the theory. 

The dependence of the magneto-phonon effect on the 
electronic Fermi energy, e^, deserves a separate discus- 
sion. It enters IIb{oj) through the filling of Landau levels 
n = 0, 1. Namely, all n+ levels with n > 2, as well as 
all levels from the upper split-off band are assumed to 
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FIG. 9. a) Raman shift and b) FWHM of the E2g phonon 
feature as a function of the magnetic field for the a~ /a^ 
(black dots) and /cf~ (gray dots) configurations. 



be empty, and all n~ levels with n > 2, as well as all 
levels from the lower split-off band are assumed to be 
filled. The n = and n = 1 Landau levels are filled 
for kz < kQi and empty at kz > /cq^i, where /cq^, kf are 
those kz at which n = and n = 1 Landau levels cross 
the Fermi level, cq = cf and ei = cf^ respectively. It 
turns out that these population effects are responsible 
for the smooth in B component of the curves in Fig. 9, 
which is different for the two polarizations (the a~ /a~^ 
curve exhibits an overall increase in energy with 5, while 
the an overall decrease). Indeed, (i) since there 

are no two 1+ and 1~ Landau levels, but a single one, 
only one of the two transitions 2~ ^ 1, 1 ^ 2+ is allowed 
for a given /c^, and the other one is blocked by the Pauli 
principle (assuming zero temperature); (ii) as the transi- 
tions 2~ ^ 1, 1 ^ 2+ contribute mainly in the <j~/<j+ 
and a~/a~^ polarization, respectively, then, depending 
on the filling of the n = 1 Landau level, the contribution 
is made in only one of the two polarizations. It turns out 
that the monotonous in B component is mainly due to 
transitions around kz = 7r/2 {H point), which becomes 
resonant with the phonon at fields near 30 T. Indeed, the 
theoretical curves in Fig. [9] are very little sensitive to the 
value of cf as it is varied from 272, the bottom of the 
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FIG. 10. a) Raman scattering spectra measured at T = 4.2 K 
(black curve) and at T = 295 K (gray curve) at B = 20 T 
in the /a^ configuration, b) Raman scattering spectra at 
selected values of the magnetic field in the cf~ /cr^ polarization 
configuration measured at T=295K. Curves are shifted for 
clarity. The vertical dashed line indicates the B=0 position. 

conduction band at /c^ = (when the n = 0, 1 levels are 
empty for all /c^), to 72 + 30 meV, but when we fully pop- 
ulate the n = 0, 1 levels at all kz by pushing k^^ ^ 7r/2, 
the tendency is reversed: the g~ j curve decreases in 
energy with 5, while the j g~ curve increases. 

D. Room Temperature Experiments 

The results presented above have been extracted from 
low temperature (4.2 K) measurements. It is however 
instructive to note that that the magneto-Raman scat- 
tering response of graphite, which in fact reflects the 
quantum character of its energy structure, survives even 
at room temperature. We illustrate this here, although 
leave a more detailed discussion of temperature effects to 
a separate report. 

Fig. [TOk ) shows the comparison of Raman scattering 
response due to i^2-,2+ transition, measured at 4.2 K 
and 295 K, in co-circular configuration. At room tem- 
perature, the L2-,2+ spectrum is weaker but still clearly 
visible. The shift of the peak towards higher energies 
and its smearing observed at room temperature is due 
to thermal population of the final state n = 2+ Lan- 
dau band. The low energy onset of the I/2-,2+ spectrum 
remains rather sharp at low temperatures since all the 



n = 2+ Landau band is empty, down to /c^ = states. 
Thermal population of n = 2+ band leads to Pauli block- 
ing of the transitions, starting from those involving the 
kz =^ final states, what quantitatively accounts for the 
shift and additional broadening of the i^2-,2+ spectrum 
at elevated temperatures. 

Remarkably, also the magneto-phonon resonance ef- 
fect in graphite survives up to the room temperature. 
Room-temperature Raman scattering spectra measured 
in cross-polarization configuration and focused on the 
phonon response are shown in Fig. [TOb). The 
magneto-oscillations of the E2g phonon are clear at room 
temperature and resemble very much those measured at 
4.2 K, presented in Fig. [51 

Rough comparison of low and room-temperature data 
points towards thermal population effects being the main 
source of the observed difference in the spectral response. 
This implies that other possible sources of spectral broad- 
ening in magneto-Raman scattering of graphite have a 
negligible temperature dependence, extraordinarily, up 
to the room temperature. 



V. CONCLUSION 

To conclude, we have used circularly polarized 
magneto-Raman scattering techniques in order to study 
purely electronic and phonon excitations in bulk 
graphite. Different series of electronic excitations aris- 
ing from i^-point carriers are observed. They are inter- 
preted using the Raman scattering selection rules other- 
wise characteristic of the bilayer graphene. The use of 
polarized Raman scattering reveals a surprisingly strong 
electron-hole asymmetry and suggests that the gener- 
ally accepted set of SWM parameters from Brandt et 
alw^ under estimates the electron-hole asymmetry. The 
magneto-phonon effect is observed in the crossed circu- 
larly polarized configuration with different amplitude of 
the effect in the two distinct crossed polarized configura- 
tions due to the contribution of the holes at the H point 
of the graphite Brillouin zone. Electronic contributions 
to the Raman scattering spectrum of bulk graphite and 
the magneto-phonon effect are observed up to room tem- 
perature. 



VI. ACKNOWLEDGEMENTS 

We acknowledge technical support from Ivan 
Breslavetz. Part of this work has been supported 
by ANR-08-JCJC-0034-01, GACR P204/10/1020, 
GRA/10/E006 (EPIGRAT), RTRA "DISPOGRAPH" 
projects and by EuroMagNET II under the EU contract 
number 228043. P.K. is financially supported by the 
EU under FP7, contract no. 221515 "MOCNA". Yu. 
L. is supported by the Russian state contract No. 
16.740.11.0146. 



11 



Appendix A: Landau levels in the four-band model 



The Hamiltonian of the A5-stacked graphite can be written as 



H 



... Vu 


V^2 


Hi 


4 


Vi\ 








... 


... 


V22 


V12 


H2 


V12 


4 





... 


... 





Vii 




Hi 






... 


... 








V22 


V12 


H2 


Vi2 


4 ••• 



(Al) 



where each element is a matrix in the Hilbert space of a single layer. Hi and H2 contain the in-plane nearest-neighbor 
coupling with the matrix element 70 as well as the diagonal shift A. V12 represents the coupling between neighboring 
layers with matrix elements 71 for the neighboring atoms, and 73,74 for the second nearest neighbors. Vu and V22 
correspond to second- nearest-layer coupling with matrix elements 72/2 and 75/2. Due to the translational invariance 
in the z direction, the problem can be reduced to that of an effective bilayer: 

Hi + Viie'^'^^ + y/ie-^^^- ^^'"■^ 
2Vi2 cos kz 



H2 + ^226 



COS kz 

2ikza- 



22^ 



-2ikzaz 



(A2) 



Here kz is the wave vector in the direction perpendicular to the layers, measured in the units of l/a^, where is the 
distance between the neighboring layers. The period of the structure is 2az, that is, two layers, so — 7r/2 < kz < 7r/2. 
Expansion of the single-layer Hamiltonian in Px^Py, the in-plane quasi- momentum components counted from the 
H — K — H line, gives 



HkAv) = 
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ri = 27iC, ^2 = 272C^ ^3,4 



273,4 
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C, r5 = 275C2 + A, C = cos kz 



a = 1.42 A is the distance between the neighboring carbon atoms in the same layer, and p± = —i%{dx i idy)- 
It is convenient to rotate the basis in the space of the 4-columns (V^i, '02, '03, '04)^ as 



^2 + ^3 
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^2 = V'l, V'3 = ^4, ^4 = 



^2 - ^3 
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so the transformed Hamiltonian becomes 
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ViP- T2 azvp+ ViP- 

ViP+ asvp- T2 -vaP+ 

V4P+ -V4P- T5-T1 _ 



(A3) 
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where V4 = v{l -\- 0^4) / \/2, V4 = v{l — 0^4) / y/2. 

In the presence of a magnetic field, described by the vector potential in the Landau gauge Ax 
the wave functions can be sought in the form 



-By, Ay = Az 



n=0 



where ^n{y) are the harmonic oscillator wave functions. When 73 is finite, the coefficient cj^^ is coupled to 

which, in turn, is coupled to Cj^'^^-' , and so on. In this situation, one can proceed either by perturbation theory or 
numerically. If we neglect 73 (set it to zero), the infinite- dimensional Hamiltonian splits into 4x4 blocks, corresponding 
to subspaces defined by 
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This greatly simplifies the calculation, as in each subspace we have a only a 4 x 4 eigenvalue problem (we denote 
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For <C Fi, out of four eigenvalues, the largest one is of the order of Fi, and the lowest one is of the order of — Fi. 
These two correspond to split-off bands and are denoted by e^++, e^--, respectively. The other two, and e^-, 
correspond to those found in the two-band low-energy approximation^. Eq. ()A9p is valid at n > 2 only. At n = 1 
there are three states corresponding to n = 1 with X3 = and energies ei++ ~ Fi — ~ — ri, and ei with low 
energy. There is one state corresponding to n = 0, with xi = = = and eo = F2. 

Optical transitions for each circular polarization are described by the operators v± — {vx ± iVy)/V^^ where Vx^Vy 
are the components of the electronic velocity operator, v = dH^^ (p)/Sp- The operator i)+ has matrix elements only 
between the n and n — 1 manifolds: 
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To obtain the polarization operator 11^ (cj), one has to sum up contributions from all kz'. 
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where each Vn contains transitions between Landau levels from the nth and (n + l)st manifolds: 
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The upper and lower ± signs in the denominators correspond to g~ j and (T+/<j~ polarizations, respectively. Both 
the electronic energies and the velocity matrix elements depend on B and k^. The Fermi wave vectors k^{B\ k^{S) 
are those k^ at which n = and n = 1 Landau levels match the Fermi energy, eo = and ei = e^, respectively. The 
sum over n in Eq. ()A12p is divergent, so we cut it off at a certain energy emax- Namely, for transitions not involving the 
split-off bands, we sum over those n for which e(^+i)+ < emax- Transitions involving the split-off bands are counted if 
^(n+i)++ < ^max- The subscqueut integration over k^ is performed numerically. The results presented in Sec. V were 
obtained with e^ax = 1 eV, and we have verified that the difference nB(co') — nB=o(^o) does not depend on emax- 
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